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We study one-flavor QCD at finite temperature and chemical potential using the 
functional renormalization group. We discuss the chiral phase transition in QCD and 
its order with its underlying mechanism in terms of quarks and gluons and analyze 
the dependence of the phase transition temperature on small quark chemical poten- 
tials. Our result for the curvature of the phase boundary at small quark chemical 
potential relies on only a single input parameter, the value of the strong coupling at 
the Z mass scale. 
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I. INTRODUCTION 

The phase boundary of Quantum Chromodynamics (QCD) is currently a very active 
frontier both theoretically and experimentally. For fixed small quark chemical potentials, 
the ground-state of QCD changes with increasing temperature from a hadronic phase with 
dynamically broken chiral symmetry to a deconfined quark-gluon plasma phase with an ef- 
fectively restored chiral symmetry. Even though the phase transition temperatures are not 
directly observable in heavy-ion collision experiments at BNL and CERN, a lower bound 
can be extracted from the experimental data These so-called chemical freeze-out tem- 
peratures can then be compared to theoretical predictions for the chiral and deconfinement 
phase-transition temperature. Since QCD is a strongly-interacting theory and long-range 
fluctuations need to be captured in order to study phase transitions, non-perturbative ap- 
proaches are indispensable for a study of the QCD phase boundary. 

On the theoretical side, various approaches are available for studies of the QCD phase 
boundary, e. g. lattice QCD simulations or functional methods. Each of these approaches 
comes with advantages and disadvantages. Lattice QCD simulations are certainly the most 
powerful tool for a study of full QCD. However, the implementation of chiral fermions 
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continues to be a non-trivial task. At finite chemical potential, the spectrum of the Dirac 
operator becomes complex, making direct lattice simulations even more difficult. In the past 
decade, however, several methods have been developed to circumvent the problems arisin^_at 
finite chemical potential, such as studies of QCD at imaginary chemicarpoteD.tial {3, [3, UQ], 

see e. g. 



Taylor expansions of the path integral or reweighting techniques jsl, Q, Isl, O, Q, 
12( 1 for short overviews. 



Refs. 
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Functional approaches to QCD, such as mean-field studies, Dyson-Schwinger Equations 
or Renormalization Group approaches, do not have problems arising from a discretized 
action or a complex-valued spectrum of the Dirac operator. However, a study of full QCD 
is not possible and a truncation of the QCD action functional is unavoidable. Therefore 
Lattice QCD simulations and continuum approaches should be considered as complementary 
approaches for studies of the QCD phase diagram. 

Dynamical chiral symmetry breaking has been studied by applying effective low-energy 
models such as the Nambu-Jona-Lasinio (NJL) model 3]. The application of these models 
is built on the assumption that QCD falls into a certain universality class, namely 0(4). 
Whether this assumption is justified or not is currently under investigation by Lattice QCD 
simulations as well as functional RG methods Q,[l5|,[l6|,0,[l8l. 

Although NJL-type models already allow to study dynamical chiral symmetry breaking 
at finite temperature and quark chemical potential, they do not contain gluonic degrees of 
freedom and they are not confining. Moreover, an ultraviolet (UV) cutoff has to be intro- 
duced in the theory. This makes the connection of these models to high-momentum scales 
and temperatures difficult. The dependence on the UV cutoff implies a parameter depen- 
dence of the model. The strategy for employing these models is usually as follows: First, 
one uses a set of parameters and the UV cutoff to fit the values of low-energy observables 
at zero temperature and zero chemical potential, e. g. to the pion mass and to the pion 
decay constant. Second, one computes the phase boundary of QCD while keeping the pa- 
rameters and the UV cutoff fixed. A shortcoming of these models is apparent: The set of 
parameters used to fit a given set of low-energy observables is not unique. Even worse, two 
sets of parameters, which both give the same results for the low-energy observables, do not 
necessarily lead to the same results for the chiral phase boundary and the location of the 
critical endpoint is not necessarily the same |l^. 

In the past few years, quite some progress has been made in connecting the low-energy 
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regime described by quark-meson dynamics with the dynamics at high temperatures, see e. 



g. Refs. 



21 



22 



23 



24 



25 



27|. Such improved models for a description of the QCD 



dynamics at finite temperature and density are mostly based on the inclusion of a Polyakov- 
Loop potential extracted from lattice QCD results. By this means, the treatment of the 
gauge-field dynamics has been outsourced while the less problematic quark-meson dynamics 
are treated self-consistently within the framework. Although all of these approaches provide 
us with a better understanding of the thermodynamics of QCD at low and high temperatures, 
they cannot get rid of the parameter dependence of the results. In addition, the back-reaction 
of the quark-dynamics on the gauge-field dynamics in terms of the Polyakov-Loop has not 
yet been fully taken into account. For a quantitative description of the QCD phase boundary, 
however, not only the gauge-field dynamics need to be taken into account: The fluctuations 
of the Goldstone modes and the radial mode beyond the mean-field approximation also play 
an important role at the phase boundary, in particular with respect to a better description 
of the susceptibilities in QCD with physical pion masses. 

In this paper, we discuss a functional Renormalization Group (RG) approach to the 
QCD phase boundary. The chiral phase boundary of the quark-meson model (bosonized 
NJL model) with two degenerate quark flavors has been studied in the local potential ap- 



proximation using a functional RG approach in Ref. [28|]. The advantage of the approach 
presented in this work is that it allows not only for dynamical chiral symmetry breaking trig- 
gered by gluodynamics but also provides access to the infrared domain of QCD dominated 
by pions. Our approach is based on ground-breaking work done by H. Gies and C. Wetterich, 



see Refs. 



29 



301 ]. There it has been shown for vanishing temperature and quark chemical 



potential that both the regime dominated by Goldstone modes and the perturbative QCD 
regime dominated by quark-gluon dynamics can be conveniently linked without fine-tuning 
using the functional RG. In Refs. 



31 



321 ] , the chiral phase boundary in the plane of temper- 



ature and number of quark flavors has been computed by studying quark-gluon dynamics 
using the functional RG. The strategy of the latter papers was to determine for which tem- 
peratures and number of quark flavors QCD remains in the chirally symmetric regime and 
thereby implicitly extracting the phase transition temperatures. In contrast, this paper aims 
to set the stage for studies of the QCD phase boundary with two and three quark flavors 
including the possibility to study the low-temperature regime and the order of the chiral 
phase transition. To this end, we use the approach discussed in Refs. 29|, l30|] and combine 
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it with the findings of Refs. 3l|, |3^: Our strategy is to follow the RG fiow starting at high 



momentum scales {p ~ Mz) down to the deep infrared regime which is dominated by the 
dynamics of Goldstone modes. This allows us to get rid of the unwanted ambiguity in the 
parameter-space as it is present in NJL-type models. Our results for the phase boundary 
will depend on only a single input parameter, namely the value of the strong coupling as at 
the initial RG scale. By this means, the scale is set unambiguously in our calculations and 
the values of all dimensionful quantities, such as the constituent quark mass or the chiral 
phase transition temperature, are eventually determined by the choice of the initial value of 
the strong coupling only. 

The paper is organized as follows: In Sect. [Tll we give a discussion of the technical details 
of our functional RG approach for a study of the QCD phase boundary. Our results for the 
chiral phase boundary at small quark chemical potentials for QCD with one quark flavor 
including a comparison to lattice QCD results are then discussed in Sect. IIIII We also discuss 



the possibility of merging our work with recent studies of t 



le deconfinement phase-transition 



in pure Yang-Mills theory using functional RG methods 33 



34| . Our concluding remarks. 



including a discussion of future extensions, are presented in Sect. HVl 



II. RG FLOW OF THE EFFECTIVE ACTION 



Throughout this paper we work in d = 4 dimensional Euclidean space and employ the 
following ansatz for the effective action for our study of the phase diagram of 1-flavor QCD: 

r = Jd'x {im + i7o/i) ^ + y [(^^)' - (^75^)1 + {d,^f + U{^') 

+ ^mr-m + + Y^iD,[A]a;rj + r„ , (1) 

where = d^6^^ — igT^^ A'^^, with being the hermitean gauge-group generators of the 
gauge group in the fundamental representation. We have introduced the shorthand {^ip) = 
ip^ipi for the color indices. In the gauge sector we have included a background gauge flxing 
term, ^ being the gauge-flxing parameter. We split up the gauge fleld into a background fleld 
and a fluctuation fleld a^, i. e. A^ = A^ + a^. The term Tgaugc contains the ghost sector 
and possible higher-order gluonic operators. We shall discuss this part of the truncation in 
more detail in Sec. Ill Fl The scalar flelds are combined in the 0(2) vector = $2) and 
we have used r = (75, i ■ 1^) in order to deflne the Yukawa interaction. The initial conditions 
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for the various couplings in Eq. ([Tj) at the ultraviolet (UV) scale A are chosen such that the 
initial effective action is given by the (classical) QCD action functional: 



fc=A 



d'x { \f;,F;, + V^i^^ 1> , (2) 



see also the discussion in Sec. IIII A[ While the inclusion of higher gluonic operators is 
discussed in Sec. Ill Fl the quark-meson part of our truncation is built around the standard 
mean-field ansatz of the effective action (large Nc ansatz; standard NJL model ansatz), i. e. 
{Z^{k = A) = 0, dtZ^ = 0; Z^{k = A) = 1, dfZ^ = 0)), which has been used extensively for 



studies of the QCD phase diagram, see e. g. Refs. |35|, l36|. Such an ansatz under 



22 



24 



25, 



ies also 



most of the recent (P)NJL studies of hot and dense QCD, see e. g. Refs. 

In the present paper, we are aiming at a dynamical connection of the high- and low- 
momentum regime of QCD. Therefore we have to go beyond the zeroth-order ansatz (stan- 
dard NJL-model ansatz) for the effective action. In the following we systematically extend 
this zeroth-order ansatz in two directions, namely in derivatives and n-point functions F*^"^ 
where n defines the number of legs. In order to study spontaneous symmetry breaking in- 
dicated by a non-trivial minimum of the order-parameter potential f/($^) in Eq. ([1]), we 
expand the potential in powers of $^ resulting in RG flow equations for the mesonic n-point 
functions, see Sec. Ill Al for details. The quality of such a systematic expansion of the ef- 



fective potential ?7($^) has been studied quantitatively in Refs. [371, l38|] and is well under 
control. On the other hand, we perform a derivative expansion which renders the n-point 
functions momentum-dependent. The latter is indispensable for a connection of the high- 
and low-momentum regime of QCD. 

Next to the zeroth-order approximation one needs to include kinetic terms for the meson 
fields in the truncation. The minimal truncation which allows for an inclusion of meson 
loops is given by the so-called Local Potential Approximation (LPA), i. e. {Z^{k = A) = 



1, dtZ^ = 0; Z^{k = A) = 1, dtZ^ = 0). This truncation has been used, e. g., in Refs. [23|, [39 1 
for a study of the quark-meson model at finite temperature and density. It indeed turns out 
that the LPA represents already a major improvement with respect to the quality of the 
critical exponents; the quality of critical exponents can be considered as a measure of how 
good the dynamics at the phase transition are captured. In the present paper, we go also 
beyond this approximation and allow for a running of the wave-function renormalizations of 
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the quark and meson fields^, i. e. (Z'0(/c = A) — > 0, dtZ^ ^ 0; Z^{k = A) = 1, dtZ^ ^ 0). Such 
a truncation renders the involved vertices momentum-dependent and improves the quahty 
of the results as it can be read off from, e. g., the quality of the critical exponents^, see e. g. 



Refs. 
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41 



42] . Aside from an extension of a given truncation with higher-order operators, 



an error estimate for a given truncation can be obtained by a variation of the regulator. By 
this means it was found in Ref. 43|] that the present truncation ([1]) gives remarkably robust 



results when applied to a study of chiral symmetry breaking at vanishing temperature and 
chemical potential. Although we have not performed such a variation of the regulator in 
the present work, it is likely that the findings in Ref. 43|] hold also in the present context of 
chiral symmetry breaking at finite temperature. This is due to the fact that chiral symmetry 
breaking sets in on scales T /k < 0.5 as we shall see below. 

We would like to point out that our truncation ([1]) is redundant since the four-fermion 
coupling Ao- is related to the scalar potential ^7($^) and the Yukawa coupling h via a 
Hubbard-Stratonovich transformation. However, this redundancy can be completely lifted 



by applying "re-bosonization" techniques 29|, l30|, |4J] which we will use here. This allows 



us to conveniently bridge the gap between quark and gluon degrees of freedom in the UV 
and mesonic degrees of freedom in the infrared (IR) regime. Moreover, re-bosonization tech- 
niques allow us to conveniently include momentum-dependent fermionic n-point functions 
(four-fermion interactions, six-fermion interactions, ^ up to arbitrary order in the RG 
flow^. In this paper, we work along the lines of Ref. |30| and give only a brief discussion of 
the "re-bosonization" -procedure in Sec. Ill El 

The effective action ([T]) has a global Ua{^) symmetry. The breaking of this global sym- 
metry is associated with topologically non-trivial gauge configurations. For the moment, 
we do not include terms that break this global symmetry, but we discuss the effect of such 
terms on our results in Sec. IIIII In one-flavor QCD, these gauge configurations play a very 



^ We would like to remark that a truncation with {Z^{k = A) ^ 0, cJtZ^ ^ 0; Z^,{k = A) = 1, dtZ^ = 0) 
represents the lowest order in the derivative expansion which allows for a dynamical connection of the 
high- and low-momentum regime of QCD, see discussion in Sec. Ill El and Refs. 



29 



3C| 



^ The derivative expansion can be continued systematically by, e. g. including terms of the form Yi;{^d^^)^ 

in our truncation, see c. g. Ref. 4C|. 
^ Note that maximal n is related to the maximal order in our expansion of the order-parameter potential 

in $^ (<f> ^ iptp) by means of the continuously performed Hubbard-Stratonovich transformations in the 

RG flow. 



7 



exposed ro 



metry 



e since they induce masslike fermion interactions which break the Ua{1) sym- 
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47 



48| . The impact of gauge-field configurations with non-trivial topology 



on the nature of the chiral phase transition of 2-fiavor QCD is not yet conclusively settled. 
In general, it is expected that such gauge-field configurations become less important with 
increasing number of quark fiavors 49|]. In this respect, we expect that dropping Ua{^) vio- 
lating terms in our ansatz for the effective action of 1-flavor QCD makes the phase structure 
more closely comparable to QCD with more than one quark flavor. 

For our derivation of the RG flow equations of the couplings, we employ the Wetterich 
equation 5ol |: 



dtT,[x] = isTr 



ri''\x]+Rk -{dA)} with T\:''>[x] 



^(1,1) I 



^ r r 1 ^ 

— rfe[x]— 



(3) 



where t = Ink/ A and A is the UV cutoff. Here, x represents a vector in fleld space and is 



defined as 



X 



X 



-q) := (Aj(-g),^^(-g),^(g),$i(-g),$2(-g)) 



and 



X = xiq) ■-- 



I AM \ 
V *2(g) J 



Thus, r[,^'^''[x] is matrix-valued in field space and so is the regulator function Rk- In this 
work, we employ a 3d optimized regulator function which is technically advantageous for 



studies at finite temperature 



51 



521 ]. The quality of such a 3d regulator in the limit 
of vanishing temperature and chemical potential has been estimated by computing critical 
exponents of 0{N) models J53|] and comparing them to those optained with an optimized 



regulator in 4d space-time [54, ISSj. Details on the regularization can be found in App. [Bl 



In the following we give the RG fiow equations in a way which does not dependent on the 
details of our 3d regularization. Reviews on and introductions to the functional RG and its 
application to gauge theories can be found in Refs. 



44 



56 



57 



58|. 
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Decomposing the inverse regularized propagator on the RHS of Eq. (jSj) into a field- 
independent and a field-dependent part, 

T^^^'\x]+Rk = Vu + Tu^ (4) 
we can expand the fiow equation in powers of the fields: 

= H"'(^/')}-H"'(MVH''(MV-- 

Here, dt denotes a formal derivative acting only on the /c-dependence of the regulator func- 
tion Rk- The powers of ^/^^k can be computed by simple matrix multiplications. The fiow 
equations for the various couplings can then be calculated by comparing the coefficients of 
the operators appearing on the RHS of Eq. with the couplings specified in our truncation. 



A. RG flow of the efl"ective potential 

In this section we discuss the effective potential. The RG fiow of the effective potential 
receives contributions from the scalar as well as the fermionic degrees of freedom: 

[/($) = [/B($) + f/ir($). (6) 

It also depends implicitly and explicitly on the gauge degrees of freedom. The implicit 
dependence affects the running of the scalar couplings whereas the explicit dependence 
would result in additional terms in Eq. ([6]). Since we are not interested in thermodynamic 
quantities such as the pressure, but only in the order parameter, we can neglect these explicit 
contributions here. The contribution of the scalar fields to the effective potential is given by 

n=—oo ^ ' 

+ Z;(^„,fe})PB(M.($))}, (7) 

where T defines the temperature and lo^ = 27mT denotes the bosonic Matsubara frequencies. 
The functions and are the wave-function renormalizations of the sigma field and 
pion field perpendicular to the heat bath. The definition of the momentum dependent 
boson propagator Pb can be found in App. \M The masses Mj of the scalar fields are in 
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general momentum dependent, = Mj(po, {Pi}), and given by the eigenvalues of the second 
derivative matrix of the potential, 

^-^•^^'^'^^ = 6i6^fi-p)) I ^'^^^^ + ^^^^(4m "'^ 5$^ = (^*r,^*D, (8) 
evaluated at the scalar background-field configuration $. In the following, we approximate 
the full potential t/ by a Taylor expansion in terms of the fields around the physical ground- 
state $0 up to quartic order^. Note that such a low-order expansion of the chiral order- 
parameter potential is incapable of describing a first-order phase transition. In particular, 
we are not able to detect the emergence a critical endpoint. However, our present work 
can be generalized straightforwardly along the lines of Refs. 61| or 28|] where first-order 
transitions have been studied within RG approaches. 

In the regime with an 0(2) symmetric ground-state ($o = 0), we use 

= lm2$2 + . (9) 

In this case, the masses of the scalar fields are given by 

A^J(^.) = 2^ + 44^|^ = ,„= + 3A,^.^ (lo) 

= 2^ M^ (11) 

and the physical masses Mj($ = 0) are degenerate. In the regime with spontaneously broken 
0(2) symmetry of the ground-state (($) = $0 7^ 0), we use the ansatz 

uu^) = ^i^' - ^ir , (12) 



which yields 



m<^) = 2^ + 4$^^^ = M^' - '^D + 2A,$^ , (13) 



= 2^ = M*^-*o)- (14) 



for the masses of the scalar fields. 



^ We neglect higher-order terms since we are not aiming at a high-accuracy determination of critical expo- 
nents, where such higher-order terms have proven their importance, see e. g. 
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The fermionic contribution Up to the effective potential U reads 



71= — OO ^ ' 

+ Zi(-j.„,{-p,})p(-)(M^)|, (15) 



where Un = {2n + l)7rT denotes the fermionic Matubara frequencies and the fermion mass 
is given by 

M| ^ Mliu^, {p,}) = ]^{h{Vn. {p^})f^^ . (16) 

The fermion propagators are defined in App. |Al As we shall see in Sec. IIIBl the wave- 
function renormalizations and Z^^ and the fermion mass have the property 

(Zi(i.„, te}))* = Z^(-z/„, fe}), {Zl{v^, fe}))* = Zli-v^, fe}) 

and (M^(i/„, fe}))* = M^(-z/„,fe}). (17) 

Thus the fermion propagator obeys 

{V^+\M^{u^,{pi}))y =V^-\M^{-u^,{pi})) (18) 

and the fermionic contribution Up to the effective potential is real-valued, as it should be: 

Uf{^) = SN^TRef^J -0^^p'{dtr^,,)Zi{u^, {p,})V^-^\M^) . (19) 

The fact that Up, and thus U, is real-valued is an important property of the effective 
potential, since it can then be expanded in powers of yU^/(7r^T^). As a consequence, we can 
expand the phase boundary in powers of /i^/(7rT)^ around fi = 0. 

The RG flow equations for the couplings m^, and the vacuum expectation value $o 
can now be calculated by projecting the RHS of Eqs. ((71) and (fTSl) onto our ansatze ([9]) 
and ( |T2l) for the potential U. In order to study the RG flow of the potential, we introduce 
the following dimensionless renormalized quantities: 

. A, 1Z^<I>1 



For the symmetric regime, we then flnd (with ^3 = g^) 

dte = (ry^ - 2)e - 8v,X^ li""^'^'^ {i, e; r/^) + 8N,v,h' ^'^^^ (t, 0,fi;vi), (21) 
dtX^ = 2r^iX^ + 20^3A5 /f ^'^'^(t, e; r/^) - SiVe^^s/i^ /f ^'('^(t, 0, jl; r/^) , (22) 
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where t = T/k and ft = fi/k denote the dimensionless temperature and dimensionless 
chemical potential. The anomalous dimensions t]^ of the scalar and 77^ of fermion field are 
given by 

vi = -d,\nZi and r/J = -a,lnZ^. (23) 

For the regime with broken 0(2) symmetry in the ground-state, we find the following flow 
equations: 

dtf^ = -{Vi + ^h + Qv, (t, 2kX^; T^i) + 2vs /f^'^'^ (t, 0; r/^) 

-8N,vSl[^^^^'\lKhMl;vi), (24) 

dA^ = 2r/^A^ + 18^;3A5/f^'^'^(t,2«:A^;r/^)+2^;3A5/f^'('^(t,0;r/^) 

-8N,v,hHi^^^^'\lKh^fi;r]^). (25) 

The threshold functions are defined in App. [B]and can be represented as Feynman diagrams 
associated with purely bosonic and fermionic loops involving the corresponding full propa- 
gators. The regulator dependence of the flow equations is absorbed into these functions. 



B. RG Flow of the Yukawa coupling 



We now turn to the calculation of the flow equation of the Yukawa coupling. Expanding 
the flow equation up to second order in the fermionic fields, we find at T = and /i = 62| : 



Of 



(2vr) 



Qo + Qo rQi + Qi 



■}) 



Qo-qo rQi- Q: 



2 



{- 



'-]) 



V^pj \m^)PbAQo + go, {Qi + gi})^ 



(26) 



where denotes the four-momenta of an incoming fermion. The flow equation for the 
Yukawa coupling h{Qo, {Qi}) is obtained from this expression by projecting it onto the 
operator -^[i/jt ■ ^ip). 

For finite temperature T, the integral in go-direction becomes a sum over Matsubara 
frequencies. Let us now discuss this integral/sum in Euclidean time direction in Eq. fl26l) . 
Since we use a 3d regulator, we can study the integral/sum in Euclidean time direction 
without discussing details of the regulator function or of the integration in spatial directions. 
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As a first approximation [62| , we set go Qo ^'iid qi — > Qi in the argument of the fermion 
mass and the Yukawa couphng on the RHS of Eq. fl26l) and take then the hmit of vanishing 
spatial external momenta Qi 0. As we shall discuss further in the next subsection, we 
also neglect a possible difference between Z^^^ and Z^j^^, thus Z^p^B = Z^^^ = We 
then obtain the following expression for 6Tj^ ^^: 

<l = / |fi9.<l({*}.««.0) + (Ps. - Ps,). (27) 

with 

^0'^^''«^^({'?^}'Qo,O)) (28) 



[MQo,0)] 



271 g2(i + ^^)2 + (^^ + i^)2 + M^iQo, 0) + r^) + {Qo - qoY + Ml ., 
and M2 = M2/Z^ and = M2/^2_ 

Let us first consider the case of finite chemical potential but zero temperature. In this 
case we observe that Eq. ( l28l) is a real number if and only if Qo is zero. This can be seen by 
writing the fermion propagator as follows: 

1 _ g2(l + r^)2 + M2(go,0) + g2-/i'-2igo/i 



+ T^Y + (go + W + M2(go, 0) \q^{l + r^)^ + (go + i^)^ + M^Qo, W 

Thus the imaginary part of the fermion propagator is linear in go and vanishes by integration 
for Qo ^ 0. This means, however, that the Yukawa coupling becomes a complex number 
for a finite external time-like momenta Qo since in this case the integrand in Eq. (l28l) is no 
longer symmetric in go. 

Now we switch on temperature. Since Qo is the Euclidean time component of an incoming 
fermion, we have Qo = (27Ti + 1)txT = Um- Taking into account that the integral in Eq. ( 128|) 
runs over fermionic momenta (go = (2n + l)7rT = z/„), we find 

Z;^Z-i(5fg^({g,},^^,0)) (30) 
°° 1 1 

Since i/^ — i/„ = 2(m — njirT is effectively a bosonic Matsubara frequency, the sum over n 
is not symmetric in n and we find that the RHS is in general a complex number for any 
given value of m, and so is the Yukawa coupling and the fermion mass. From Eq. ( 130|) . we 
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conclude 

(/i(z/„,0)r = /i(-z/^,0), (31) 

and equivalently for the fermion mass. Thus we have found that the Yukawa couphng for 
a given external momenta is in general complex-valued when evaluated at a finite quark 
chemical potential. This is not an issue as long as we take the full momentum dependence 
of the Yukawa coupling into account, as we have argued in Sec. Ill A[ In the following we 
shall restrict ourselves to a momentum-independent Yukawa coupling which has been suc- 
cessfully employed in studies of the quark-meson model with two quark flavors at vanishing 



chemical potential, see e. g. Ref. 



41 



42 



621 ]. This requires care in finding a proper ap- 



proximation scheme that gives us a real-valued effective potential U without computing the 
full momentum dependence of the couplings. The idea for constructing such a scheme is to 
expand the theory around the limit vrT/A; ^ of vanishing external momenta. This might 
appear dangerous for k < ttT, but what comes to our rescue is the fact that TcT/k < 1/2 
above the scale k^sB at which QCD enters the chirally broken regime. For scales k < k^^sB, 
the fermions acquire a mass due to the presence of a quark condensate and therefore the 
fermionic contributions to the flow decouple rapidly anyway. Thus, once chiral symmetry 
is broken, our approximation in the fermionic subsector should not influence our results 
much. For the purpose of implementing this truncation scheme, we introduce the following 
dimensionless quantities: 



X 



""""T ^^,B = -^- (32) 



Now we rewrite Eq. (130|) in terms of dimensionless propagators: 

Z7Z^i(A:5fg^(te},z>„,0)) (33) 

1 1 



^^^'''"'°^^'\5^x2(l + r^)H(z>„,+i/2)Hm2(z>„,0)a;2(l + rB) + (i>™-z>„,)Hm|^ 
Assuming that Um is a small parameter, we can expand the boson propagator in powers 

of Um'. 

1 



X2(l + Tb) + [Vni - ^nY + 



B 

x\l + rB) + ul + ml\ x\l + rB) + K + 'ml 
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Equivalently, we expand the Yukawa coupling and the masses: 

h{Dm, 0) = h + h^^^vm + 0{vl) , (35) 
'm^,B{i>m, 0) = m^^B + rn-l^B^m + C'(t'm) • (36) 
Keeping only the zeroth order in these expansions, we obtain for 5T^ ^^{{^i]^ ^m-, 0): 

(fc5fg^({g,},z>„,0) 

oo 

= /^s^ \^ t t. f37) 

^ a;2(l + r^)2+ (^„ + i/i)2 + m2 2;2(^1 +^^) + j>2 + ^|^^ ■ v i 



n=— oo 



We observe that this expression is a real number for all fx. Thus a Taylor expansion of this 
expression around fl = generates only terms with even powers in fl. 

Inserting Eq. fl371) into Eq. fl27j) and incorporating the gluons in the same way as discussed 
here for the scalar fields, we obtain the final result for the flow of the Yukawa coupling: 



-32^3//i'C2(iV.){LS5^)'(')(t, Kh', fi, 0; r/^, ry^) -l^£S^^^)'(^)(t, f,h', /i, 0; r/^, vf)} , (38) 

where rj^ = —dtlnZp. Moreover, we have {ml = e,ml = 2n\^) in broken regime and 
ml = ml = e in the symmetric regime. The threshold functions associated with triangle 
diagrams are defined in App. [Bl We have checked that our results agree with those in 



Ref. 30| in the limit T — ^ and — if we use a four- dimensional regulator function. 



C. RG Flow of the wave-function renormalizations 

The flow equations for the fermionic wave-function renormalization can be extracted from 
the RHS of Eq. fl26|) along the lines of the calculation of the Yukawa coupling. Projecting 
Eq. ( !26l) onto ilj{—'jiQi)ip and then taking the hmit Qo ^ and Qi 0, we obtain the flow 



4t;a 



= ^C2(7Vc)/|4Al£''^'^'^(t, Kh', /i, 0, 0) - 8(1 - /i, 0, 0) 

+f (1 - (A/'i^r'^'^(t, f,, 0, 0) + M!:f'^'\i, Kh\ f,, 0, 0) 
+ ^h4 (^™(t;/./.^/i,m^,0) + ^^™(^)(^;«:/^^;i,m^,0)) V (39) 
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where the corresponding threshold functions can be found in App. [Bl Note that we do not 
display the dependence of the threshold functions on the anomalous dimensions t]^ of the 
scalars, 77^ of the fermions and rjF oi the gluons for brevity, but we take it into account in 
the numerical evaluation of the flow equations. We find agreement with the equation for rj^ 
in Ref. 30|] in the limit T — > and /i — * if we use a four-dimensional regulator function. 
A flow equation for z|),((5o! {Qi}) could in principle be obtained in the same manner by 
projecting Eq. fl26l) onto i'i—'JoQo)'^- 

The derivation of the scalar wave-function renormalization can be performed along the 



lines of Refs. 



30 



37 



62l |. Projecting the flow equation onto and taking the limit of 



Qo = and Qi = for the external momenta, we find the wave-function renormalization Zi 



+ ^-^NXnMf^^^'\lKh\n:4)- (40) 

The threshold functions are defined in App. [Bl For vanishing temperature and quark chem- 
ical potential, we find that the equation for r/^ agrees with the equation for 77^ provided we 
employ a four-dimensional regulator function. 

Here and in the following we neglect that Z^^ 7^ Z^ and 7^ Z^ and work in the 
approximation Z^^ = Z^ and Z^j^ = Z^. For our purposes, this is justified since we are only 
interested in chiral symmetry breaking in QCD but not in a calculation of thermodynamical 
quantities above Tc- For scales k > k^sB, "we have zf ^ Z^ since T/k < 1. For scales 
/c -C T, we approach the three-dimensional limit and the RG flow is driven only by the 
lowest Matsubara modes. Since the lowest Matsubara frequency for the bosons is zero, 
the dependence of the propagators on z|| drops out, see App. |X] for our definition of the 
propagators. In case of the fermions the lowest Matsubara frequency is proportional to the 
temperature T. Therefore the fermions effectively decouple from the RG flow for <^ T 
and our approximation Z^ hardly affects the RG flow. Overall, the distinction of 

zl and Z^- plays only a quantitatively important role for temperatures T > T^. There the 
mid-momentum regime is not protected by a dynamically generated mass gap from chiral 
symmetry breaking and modes with k ^ T can actually probe the difference between zf 
and Z^. 
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D. RG Flow of the four-fermion interaction 

The flow equation for tfie four-fermion interaction can be obtained as well by projecting 
the expansion ([5]) of the RG flow equation onto our ansatz ([T]) for the effective action. In 
anticipation of what follows, we note that we only need to take contributions arising from 
the fourth order term in the expansion ([5]) into account. These are contributions from so- 
called one-particle irreducible (IPI) "box" -diagrams. As we shall see in Sec. Ill Et we do 
not need to compute IPI four-fermion self-interaction diagrams (~ A^) and so-called 1 PI 
"triangle" -diagrams (~ Xaff' and ~ \ah?), even though these diagrams would contribute to 



the RG flows of the four-fermion couplings in a non-rebosonized study 3ll, l32, |43|, l63|] . As 
a consequence, it is sufficient to consider the limit Ao- — > on the RHS of the flow equation 
of A^. We find: 



-^K^'+^ry) (41) 



with 



Pt = ^^^i^L\^^l\^'\lnh\j.r^,^^^^^ 1, (42) 

PC = -'-^^^^f -h\ ;2, /i, 0, 0) + LS-;f (t, .h\ 0, 0) j . (43) 

The threshold functions can be found in App. [Bl We do not display the dependence of the 
threshold functions on the anomalous dimensions of the corresponding fields for brevity but 
we take it into account in the numerical evaluation of the flow equations. 



We have chosen the same Fierz transformations in the Dirac algebra as in Refs. [29|, |30 |. 
In the present study we discard additional four-fermion interactions of the type [ip'joip)^ 
which are generated in the finite-temperature RG flows. However, such interactions are 
suppressed for scales k > T compared to the included four-fermion interaction anyway. We 
have also neglected four-fermion interactions, such as a vector-channel, in our truncation of 
the effective action. Therefore our results for the phase boundary will depend slightly on our 
choice of Fierz-transformation with respect to Dirac and color indices. However, it has been 
checked in Ref. js^ that results for low-energy observables at zero temperature obtained in 
different Fierz decompositions involving a color-singlet scalar-pseudoscalar channel agree on 
the 1% percent level. We would like to stress that it is possible to fully resolve such a Fierz 
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ambiguity in larger truncations within the 



"re-bosonization" techniques are apphed 



29 



unctional RG approch 3l|, |32, |43|] even when 



63 



64| 



E. RG Flow of the rebosonized couplings 

Let us now briefly discuss the so-called "re-bosonization" procedure which we apply in 
order to resolve the redundancy in our ansatz for the effective action ([T]). The redundancy 
originates from the fact that a Yukawa coupling together with a bosonic potential can be 
transformed into a four-fermion interaction and vice versa. In order to lift this redundancy 
we work along the lines of Ref. |30| and allow for /c-dependent scalar fields and $2,fc- 
The flow equation ([3]) changes then as follows 0, Q : 





f d\ j 


+ 


J (2vr)4 I 



AF \ 
•9**1,^ + T^^**2,fc , (44) 



where the first term on the RHS is simply the flow equation ([3]) for fixed fields and $2,fc 
and the second term takes care of the fact that the scalar fields change under a variation of 
the RG scale k. In the following, we span the RG flow of the scalar fields and $2,fc by 
the corresponding field itself and a fermionic composite operator with the same quantum 
numbers: 

at$i,fc = (^757/;) dtak + ^i,kdtPk , (45) 

dt^2,k = ^ i^^p) dtak + ^2,kdtf3k . (46) 

The functions ak and Pk determine the transformation of the scalar fields under the RG 
flow and can be derived unambiguously from enforcing several conditions: The flow of 
Ao-(q'^) must vanish on all scales k and for all g^, the Yukawa coupling must be momentum 
independent and the flow of the wave-function renormalization of the scalar fields must obey 
dtZ^{q'^ = k"^) = —rj^Z^, see Refs. Q for details. Using the initial condition A^-lfc^A = 
for the four-fermion coupling at the UV cutoff scale A, the first condition ensures that no 
coupling Ao- is generated in the RG flow. We stress that it is this "re-bosonization" procedure 
which allows us to bridge the gap between the perturbative quark-gluon regime in the UV 
and the regime dominated by massless Goldstone modes in the IR without performing any 
additional fine tuning. From a technical point of view, this technique allows us to include 
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(momentum-dependent) four-fermion interactions up to arbitrary order, provided we do not 
truncate our ansatz for the scalar potential f/($^). 

By applying the field transformations (H5l) and (H6l) . we modify the fiow equations for 
scalar couplings. We find for the fiow equations in the symmetric regime: 

+ 1 



dte = dt.e 



((i+e)g. + i) [piy+p^y) 



+ 2 ((1 + efQ, + 1 + 2e) (^^f/ + P^h 
+ 4^ (1 + e) (1 + (1 + e)g.) + f3lh 



(47) 
(48) 

(49) 



Similarity we obtain the following set of fiow equations in the regime with broken 0(2) 
symmetry: 



dfK = dfK 



((l + «:A^)g. + l)(/?f/ + /3£/i^ 



+ 2 ((1 - KX^fQ, + 1 - 2kX^) [[3f/ + /?f>^) , 

+ 4^ (1 - kx^) (1 + (1 - kx^)q.) + Plh^) 



(50) 
(51) 
(52) 



The function occuring in the equations for the symmetric and the broken regime measures 
the suppression of the four-fermion interaction for large momenta, which we treat in an s- 
channel approximation. It is defined as 13 • 

dtCK{k^ T, /X, M^) - A^(0, T, fi, M^)) 



(53) 



atA,(0, T,/i,M^) 

Note that Q^, depends on the temperature T, the quark chemical potential /i and the quark 
mass M^. In order to compute Qa-, we would in principle need to compute the full mo- 
mentum dependence of the four-fermion interaction. For simplicity, we do not perform an 
explicit computation of the momentum dependence but model it with the aid of theoretical 
constraints. First, we assume that < in order to be consistent with unitarity at T = 0. 
Second, the four-fermion interaction in the s channel can be considered to be roughly point- 
like once the quarks acquire a mass. At finite temperature, the quarks acquire an additional 
thermal mass which further suppresses the momentum dependence. According to Ref. 30| . 
we therefore model Qa by employing a threshold function that captures these constraints: 



(54) 
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Here, Q° is a negative constant at our disposal. We choose Q^^ = —1, but we have checked 
that our resuhs for the phase boundary and in particular for the curvature at small chemical 
potential change only on the 1% level when we vary from (5° = —1 up to = —0.01. 



F. Running of the gauge coupling and gluonic anomalous dimension 



Finally we need to discuss the running of the gauge coupling which is one of the key 
ingredients of our study of the QCD phase boundary. From now on we restrict our discussion 
to Landau-gauge QCD. 

In this work, we use two ansatze for the running of the coupling in Landau-gauge QCD. 



This allows us to give a theoretical error estimate for our results. T 
extensively discussed at both zero and finite temperature in Refs. 



le first ansatz has been 
3^ 



31 



65| . It is based 



on the following truncation in the pure gluonic part of the effective action: 



(55) 



Such a calculation of the coupling employs 



RG framework 



44 



65 



67 



68 



69 



7G 



71 



the background-field formalism 66| within the 



72 



73|. 



The ansatz (15511 for the pure gluonic part of the effective action used for the determination 
of the running coupling includes an infinite power series of the gauge-invariant operator 
F^v^'^u- The truncation includes arbitrarily high gluonic correlators projected onto their 
small-momentum limit and onto the particular color and Lorentz structure arising from 
powers of It represents a gradient expansion in the field strength to arbitrary order 

but neglects higher- derivative terms and more complicated color and Lorentz structures. 
Using the background-field method, the /^-function of the running coupling g is related to 



the wave-function renormalization of the background field 



66| via 



C^tttFE = '7FE "EE with 



'7FE 



-ainz 



(1) 



EE 



(56) 



which is a consequence of the non-renormalization of the product of the background field 
and the bare coupling. The coefficient of the first term Z^^ = Zp^/A in the effective action 
fl55|) evolves with the renormalization scale k and is successively driven by all other operators 



in the action. In Refs. 



31 



65|, the authors keep track of all contributions from the flows 



of the Z^*^ to the flow of the running coupling. An infrared fixed point for the running 
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coupling at zero temperature has been found in Ref. [65| with dp^iT = 0) G [5.7,9.7]. The 
uncertainty arises from an unresolved color structure in the calculation. In the following we 
do not compute the running of the coupling from the truncation (l55l) explicitly but use the 



results from Refs. 



31 



32j with a*pp{T = 0) = 5.7. 



One of the main findings in Refs. 3ll, l32] is that the coupling exhibits a non-trivial 



infrared fixed point at finite temperature which has been recently confirmed by Lattice 
QCD simulations [3]. In the low momentum regime, the solution of the RG equations 

ope determined by the infrared fixed point a^^ of the 



exhibits a linear behavior with a s 



spatial 3d Yang-Mills theory [3ll, l32| : 



apE{k < T) 




(57) 



The value of the infrared fixed point is given by a^^ ~ 2.7. The actual presence of this 
finite infrared fixed point is important for temperatures around the chiral phase transition 
while the actual value of a^^ is of less importance for a study of the chiral phase transition 
temperature^. Indeed, the running coupling obtained from the truncation ( 155|) has been 
successfully used to determine the chiral phase boundary of QCD in the plane of temperature 
and number of massless quark flavors 3l|, 1321]. 

Since our work relies partly on the background-field method, we would like to dis- 
cuss briefly its advantages and disadvantages. The application of the background-field 



method to functional RG flow equations has been proposed in 
in Refs. 



44 



65 



67 



68 



69 



70 



71^, 



72 



671 ] and further developed 



731 ]. The background-field method provides a conve- 



nient framework for a study of gauge theories since it allows us in principle to construct 
a gauge-invariant effective action in a comparatively simple manner. The background- 
gauge fixing procedure 66| together with the regularization lead to regulator-modified Ward- 



Takahashi identities (mWTI) 



68 



|69|, 170|. Here, we only employ an approximate solution to 
31 . 3^ . To be more specific, we identify the 



the flow in the gauge sector as obtained in Refs. 
RG flows of the background field with those of the fluctuation field; for a treatment of the 
difference of both we refer to Ref. 7l|. The identification of the background and fluctuation 



field results in a flow which is no longer closed and which does not satisfy all constraints 



The actual value of may play an important role for the study of bulk thermodynamic quantities such 
as the pressure at high temperatures. 
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from the mWTI [4^, |67|, l68|, l69|, 170|, lDJ, 172|, |73|. In this work, we assume that the loss of 
information due to the identification of the background and the fluctuation field as well as 
corrections due to the mWTI are quantitatively small in the region of physical interest and 
do not severely affect our results, see Ref. 72] • The advantage of our approximations is 
that we obtain a gauge-invariant approximate solution of the theory. In the following we 
work in Landau-deWitt gauge whenever the background- field method is involved^. Since 
we shall also employ the coupling from lattice QCD in Landau gauge, our analysis of the 
phase boundary provides some quantitative insight into the quality of our approximations 
involved in the gauge-sector when treated within the background-field formalism. 

In order to "measure" the impact of the gauge field dynamics on the QCD phase boundary 
and to estimate the theoretical error of our results arising from the truncation in the gauge 
sector, we also employ the running of the gauge coupling in Landau-gauge as measured 

771 ]. In Landau-gauge QCD, the running of the gauge coupling at 



on the Lattice in Ref. 



vanishing temperature has been cornputed using lattice simulations [77 



Dyson-Schwinger equations 82 



83 



85l ] and functional RG methods 



78 

86, 



80 



8l|, 



79. 



87|. It can 



be defined by means of the ghost and the gluon propagator which is a consequence of the 



non-renormalization property of the ghost-gluon vertex 



82 



83 



88j: 



(58) 



47rZA(T = 0,p2)z2(T = 0,p2)' 

where Za,c denotes the dressing functions of the gluon and the ghost, respectively. The 
momentum dependence of the dressing functions are characterized by a power-law behavior 



in the deep IR 



-2Kr 



2\Kc 



(59) 



The exponents are related by the Landau-gauge sum rule in ci = 4 dimensions [89|, |90|, |9l| . In 
this work, we have suitably amended the lattice propagators in Ref. 77| by their perturbative 
behavior in the ultraviolet regime and the corresponding power laws in the IR. This yields 
an infrared fixed point as{T = 0) ~ 2.3. 

For our finite temperature studies, we have adapted the running of the gauge coupling 
f l58|) such that it is governed by an infrared fixed point for momenta p < 27rT according to 



^ Note that Landau gauge is known to be a fixed point of the RG flow 



75, 



76| 
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the results in Refs. [3l|, |32| |: 



«Rcf.(P « T) ^al,^ + ai (I) + a2 (I) + • • • • (60) 

Here we drop all higher terms and choose a*^^ = 1 and determine Oi and a2 such that 
the coupling (158!) and its derivative with respect to p are connected continuously with the 
ansatz fl60l) at the scale set by the lowest non- vanishing bosonic Matsubara-mode uj- = 2'kT. 
Although the actual values for ai, 02 and at^^ may differ from the values chosen here, the 
arising uncertainties for the QCD phase boundary can be estimated by a comparison with 
the results obtained from the coupling defined in Eq. (l56l) . In any case, the question whether 
the ground-state of QCD is governed by chiral symmetry breaking or not is controlled by 
the running of the coupling in the mid-momentum regime (0.5 GeV ^ p ^ 1.5 GeV) as we 
shall see below. In this momentum regime, we have aRef. > ctpE- 

From now on we identify the scale k"^ set by the cutoff function with the momentum 
scale p^. This nontrivial assumption is justied, because the regulator function which enters in 
the calculation of the running coupling specifies the Wilsonian momentum-shell integration 
in such a way that the RG flow of the coupling is dominated by fluctuations with momenta 
p^ k"^ . For our calculation of the phase boundary, we need not only the strong coupling, 
but also the anomalous dimension rjF which we estimate from 

dtOtKcL {p^ = k^) = r/Rof. ttRcf. (p^ = /c^) . (61) 

By estimating the gauge-field contributions to rip = —dfZp from ?7Ref. and by using ORef. in 
our calculations, we assume that the running of the coupling as found in Landau-gauge QCD 
can be identified with the running of the coupling found in Landau-DeWitt gauge within the 
background-field formalism. This means that we ne glec t possible differences between the 



RG flows of the fluctuation and the background field [71|, [72] • The two-loop running of the 
coupling is indeed identical for both gauges despite the fact that the couplings are defined 
differently. In the deep IR, there is qualitative agreement between the Landau gauge and the 
Landau-DeWitt gauge indicated by the presence of a non-trivial attractive IR fixed point. 



This suggests a deeper connection between both gauges [33|] and justifies the use of both 
results for the coupling for our computation of the QCD phase boundary. The differences in 
the results can then be considered as a measure of the influence of the gauge-field dynamics 
on the phase boundary. 
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Finally, we need to discuss how the quarks affect the running of the strong coupling. In 
this work, we only use a one-loop RG improved quark contribution to the gluonic anomalous 
dimension and therewith to the running of the strong coupling. This contribution can 
be straightforwardly computed from the effective action ([T]). Using the regulator shape 
functions (1B4I) . we obtain'' 

,~ Nf ( 1 1 \ 4 ^2 

where denotes the dimensionless quark mass, and t and fi denote the dimensionless 
temperature and quark chemical potential, respectively. We recover the standard one-loop 
result in the limit of vanishing quark mass, chemical potential and temperature, which is as 
it should be. The gluonic anomalous dimension rjA is given by simply adding the gluonic 
contribution and the quark contribution: 

Rcf./FE , 
Vf = VRef./FE + Vq ■ (63) 

The running coupling of the strong coupling g including the back-reactions of the quarks is 
then obtained by solving the differential equation 

a Ref./FE /^.N 

C't«Ref./FE = Vf «Ref./FE , (64) 

which is coupled to the the RG flow equations for the Yukawa coupling h and the minimum 
of the scalar potential n. 



III. THE PHASE BOUNDARY OF 1-FLAVOR QCD AT SMALL CHEMICAL 

POTENTIALS 

A. Initial conditions and fixed-point structure of QCD 

Let us now discuss our results for the phase boundary of QCD with one quark flavor. As 
initial conditions we use the running coupling as measured at the Z-boson mass scale Mz, 



as{Mz) ~ 0.118 92|, and set the renormalized Yukawa coupling to h{Mz) = 0.01, the 



^ The quark contribution rj^ to the gluonic anomalous dimension should not be confused with the anomalous 
dimension associated with the quark propagator. 
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renormalized scalar mass to m?{Mz) ~ 5M| and the renormalized four-boson coupling to 
= 0. With this choice, the effective action at the initial scale Mz is effectively given by 



k=Mz 



j rf^x|lF;,F;, + ^i^^| . (65) 



Our choice of initial conditions is such that the results for the phase transition temperature 
depend only on the value of the strong coupling at the initial RG scale. The choice of the 
value for strong coupling at the initial scale eventually determines the absolute values of 
observables such as the phase transition temperature or the constituent-quark mass. Let 



301 1 and then generalize them to the case of 



us briefly recapitulate the results from Ref. 
finite-temperature QCD. 

The universal features of spontaneous chiral symmetry breaking in QCD can be nicely 
illustrated in terms of the fixed-point structure of the scalar couplings. This structure can 
be considered as the analogue of the fixed-point structure of four-fermion interactions in a 



purely fermionic language (Sll, l32]- For this purpose, we define the coupling 



The definition of this coupling is simply motivated by the relation Ao- ~ h'^/m^ between the 
four-fermion coupling X„ in the gauged NJL model and the scalar couplings rh^ and h"^ in a 
bosonized gauged NJL model. 

The flow equation for the coupling e can be straightforwardly derived from the flow of 
the coupling e and the Yukawa coupling h. We find* 

-(2- ^2v,C,{K)g' |LS;;^)'(^)(t,0,/i,0) - l4r)'(^)(t,0,/i,0)}) I 

-2(piy+p^y)i^. (67) 



We have neglected the anomalous dimensions for simplicity since here we are interested only 
in the weak coupling regime of QCD. We find that l3i has a UV repulsive fixed point e* and 
an IR attractive fixed point 62 if the gauge coupling g is smaller than a critical value Qct, see 
Fig. [TJ For g = Qcr, the fixed points annihilate and the flow of e is not bound by any fixed 



Here we drop all arguments of the threshold functions associated with anomalous dimensions. 
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Figure 1: Sketch of the function of the couphng e. The figure illustrates the dependence of 
the flow of e on the gauge coupling g. The /?- function has an UV repulsive fixed-point e| and IR 
attractive fixed-point £2 if the gauge coupling is smaller than a critical value g^j.. The fixed 
points annihilate for g^ = g1^ and no fixed points are present for g^ > g^^. The red dashed line 
illustrates the effect of finite temperature on the result for g"^ < g1^ indicated by the black line. 

points for g > Qcr- In the latter case, the system flows towards the regime with broken chiral 
symmetry characterized by a negative scalar mass parameter e = m'^/k^. Thus, g > Qct 
represents a necessary condition for chiral symmetry breaking and the question of whether 
the QCD ground-state is chirally symmetry or not has been traced back to the strength 
of the gauge coupling g relative to its critical value gcr- We would like to mention that 
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431] where chiral symmetry 



these considerations fully correspond to those in Refs. 
breaking has been studied in terms of four-fermion interactions. The critical value of gcr can 
be computed analytically at T = and /i = in the limit e ^ 1 from Eq. ( 1671) : 

"- = 7^^7^tI?T ^ iV,(72^ - const, for iV, > 1 . (68) 

Air C2{J^c) 

The numerical factor arises from the evaluation of the threshold functions^. Since the influ- 
ence of the Yukawa coupling and the four-boson coupling can no longer be neglected near 
the chiral transition transition scale^°, this can only serve as an estimate. 



^ The result for acr deviates from the value in Ref. [3Qi| due to the choice of a 3d optimized regulator- 
function instead of a 4d optimized regulator-function. The difference between both results is expected to 
be smaller in the case of finite temperature since the 4d and 3d optimized regulator function coincide in 
the limit T/fc 00. 

^° With respect to the influence of the scalar couplings, the estimate for g^r given here does not necessarily 
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The fixed-point structure of the couphng e allows us also to divide the initial conditions 



into different sets. As in Ref. |30|, the QCD starting point is given by eA = ^ ^ 1 since 
the scalar field is then an auxiliary field at the initial scale A = Mz and its wave-function 
renormalization is tiny. Our choice for the initial conditions corresponds to a value of e 
which is close to e^- Starting the RG fiow from a set of initial conditions obeying e > 63 
at the UV scale A = M^, the system fiows into the fixed point 63 which then controls the 
evolution over a wide range of scales. This explains that the IR physics at zero temperature 



3- 



are not sensitive to the choice of the initial conditions 

At finite temperature, the situation changes slightly. For a given value of the gauge 
coupling the depth of the minimum of the /5e-function is increasing with an increase in 
the dimensionless temperature T/k and the distance between the fixed points increases, see 
Fig. [H We find 

lim 1*2 00 . (69) 

t— >oo 

Thus the critical value for the gauge coupling g increases with increasing T/k as well. 
This can be understood phenomenologically: The formation of a quark condensate requires 



stronger interactions since the quarks are thermally excited [3l|, l32[ . In order to leave the 
result for the chiral phase boundary unaffected by the choice of the initial conditions, we have 
to choose the initial scale A in such a way that the flow is still initially governed by the flxed- 
point £2 given at T = 0. In practice, this can be achieved by choosing a large value for the 
initial scale A such that T/A <^ and /i/A <^ for the temperatures T and quark chemical 
potentials /i under consideration. Our choice A = Mz translates into T/A O{10~^) and 
/i/A ~ 0(10^'^) for all of our numerical evaluations of the flow equations. Thus, our results 
for the phase boundary are not contaminated by the choice of the initial conditions. 

Our discussion shows that the scale for all dimensionful quantities is set by the interplay 
between the perturbative running of the gauge coupling with its initial value, the critical 
value gcT, and the existence of the flxed point i^, provided we chose QCD-like initial con- 
ditions. The presence of the flxed point €2 ensures that the actual initial values for the 



even if one neglects that different regulator 



have to agree with the value found in Refs 
functions have been used there. In Refs. 31,|32, |43[, the authors study the influence of gluodynamics on 
the quark dynamics in QCD without using rebosonization techniques. The flow of the scalar couplings 
is then completely encoded in the RG flow of the four-fermion couplings and the estimate for within 
such a "pure" quark-gluon framework should be considered to be more accurate. 
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scalar couplings do not affect the absolute values for the constituent quark mass or the 
phase transition temperature, for example. The initial value for the gauge coupling g and 
its (logarithmic) running, in combination with the critical value QcriT/k, fi), then set a scale 
kcr(T,fi). For given values of T and fi, it is given by the solution of the equation 

9ih f ) = f ) . (70) 

The scale kcr is naturally related to the scale Aqcd at which the gauge coupling becomes 



large. Thus the scale for all dimensionful quantities is se ; 
conditions at the UV scale. Along the lines of Ref. 31 



3y Aqcd, independent of the initial 



32l ]. one can actually use Eq. ([70 



to determine an upper bound for the chiral phase boundary. This is done by seeking the 
lowest temperature for a given /i above which Eq. (!70|) does not have a solution anymore. 
We shall discuss this approach further in the next subsection. 

Let us finally compare our approach with (P)NJL-type models. In the case of (P)NJL- 
type models, it is necessary to introduce a UV cutoff A which is usually on the order of 
1 GeV. At this UV cutoff scale the gauge degrees of freedom are considered to be integrated 
out. The cutoff can therefore be considered to have a physical meaning and is needed to 
define the theory. The strategy is then to choose the value for the four-fermion couplings 
(or the scalar couplings in the context of the bosonized NJL-model, i. e. the quark-meson 
model) which reproduce the values of low-energy observables, such as the pion mass, the 
constituent-quark mass and the pion decay constant. These values of the couplings are then 
used to compute the chiral phase boundary in the plane of temperature and quark chemical 
potential. The shortcoming of this procedure is that it does not result in unique values 
for the initial values of the couplings at T = /i = 0: Different sets of initial values for the 
couplings can reproduce the same values for the low-energy constants under consideration 
but lead to different predictions for the chiral phase boundary and the location of the critical 



point^^, see e. g. Ref. 19|. If we interpret (P)NJL-type models as a low-energy formulation 
of QCD, we can understand their parameter ambiguity in terms of the fixed point structure 
of the coupling e by considering the limit of vanishing strong coupling g. In this limit, the 



11 



Even though the predictions from Lattice QCD simulations for the curvature of the phase boundary at 
small chemical potentials and the location of the critical endpoint do not yet agree as we 
would like to stress that the reasons for these differences are completely different from those encountered 
in the context of (P)NJL-type models. 
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fixed-point value corresponds to the value of the inverse of the critical coupling in the 



NJL model 
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301 ] ■ Choosing initial conditions with eA < e*, the system is driven by strong 



four-fermion interactions and flows towards the regime with broken chiral symmetry. In 
our approach, we do not encounter such an ambiguity in the choice of the initial conditions 
because of the dynamically included gauge degrees of freedom and the presence of the fixed 
point €2- The scale for the IR physics is uniquely fixed by our choice for the initial value of 
the gauge coupling g and all absolute values of the observables should be interpreted in the 
light of this scale-fixing procedure. 

In the context of (P)NJL-type models, one might be concerned that one needs to adjust 
the initial values of the couplings at finite temperature. This concern is based on the 
observation that the temperature effects at the UV scale Anjl might be non-negligible since 
^/Anjl ^ 0.3 for typical UV-cutoff scales Anjl and temperatures used in (P)NJL-type 
models. As we have argued above, this problem is not present in our approach. 

Despite the nice features of our RG approach, the present truncation is by no means 
complete, in particular with respect to additional operators which affect the mid-momentum 
regime where the transition to the regime with broken chiral symmetry takes place. In this 
regime, the flow is sensitive to higher-order operators owing to strong coupling. Examples 
we have in mind here are the inclusion of the Polyakov-Loop in the way proposed in a 
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341 ]. or the inclusion of operators 



RG framework for pure Yang-Mills theory in Refs. 
associated with instanton effects. The latter are expected to be particularly important for 
1-fiavor QCD. The systematic errors arising from neglecting such operators are discussed in 
the next subsection. In any case, we think that our approach is promising and sets the stage 
for future works in this direction. 



B. Results for the curvature of the phase boundary of 1-flavor QCD 

Let us now discuss our numerical results for the (chiral) phase boundary of QCD with 
one quark flavor. As a first result, in Fig. [2] we show the temperature dependence of the 
quark mass for vanishing quark chemical potential. The quark mass is related to the order 
parameter n by = nh"^ and tends to zero continuously for increasing temperature, which 
indicates a second order phase transition. This is expected since so far we have neglected 
Ua{^) symmetry breaking terms in our study. The results shown are obtained from the 
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Figure 2: Temperature dependence of the quark mass for vanishing chemical potential as obtained 
from a calculation employing ctRef.- The quark mass goes continuously to zero at the critical 
temperature indicating a second order phase transition. 

calculation employing the strong coupling a-Rci.- For the quark mass at T = and /i = 0, 
we find 

mJ"^-(T = 0, /i = 0) = 434 MeV and m™(T = 0, /i = 0) = 430 MeV (71) 

for aRcf. and apE, respectively. For T = and /i = 0, the corresponding scale kcv at which 
the RG flow enters the regime with broken chiral symmetry is given by 

A;^/- (T = 0, /i = 0) = 470 MeV and A;™(T = 0, = 0) = 330 MeV. (72) 

For the phase transition temperature, we find 

Tj^«f-(/i = 0) = 110 MeV and T™(/i = 0) = 76 MeV (73) 

for ORcf. and apE? respectively. The difference in the phase transition temperature is 
mostly due to the differences in the running of the gauge coupling in the mid-momentum 
regime (0.5 GeV ^ p ^ 1.5 GeV). Note that the RG flow at flnite temperature is less sensi- 
tive to the signiflcant difference between the IR flxed point values of the couplings aRef. and 
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ttpE at vanishing temperature^^. This is because the RG flow of the couphng for k < 2ttT 
is governed by the fixed point of the underlying 3d Yang-Mills theory. Moreover, chiral 
symmetry breaking is triggered in the mid-momentum regime rather than in the deep IR 
regime. Therefore we would like to add that we expect in general that our results do not 
strongly depend on the behavior of the coupling in the deep IR, independent of the fact that 
the IR running of the strong coupling at finite temperature is governed by the underlying 3d 
Yang-Mills theory. To be more specific, the gluons strongly influence the matter sector in 
the mid-momentum regime via gluon-induced four-fermion interactions. As a consequence, 
the gauge degrees of freedom drive the quark sector to criticality depending on the actual 
temperature, see Subsec. IIII Al for details. Below the chiral symmetry breaking scale fccr, the 
quarks acquire a finite mass. In the deep IR (fc ~ p < 200 MeV), where currently debated 
differences between a scaling and a decoupling solution in the gauge sector become signifi- 
cant 93| , the quarks are decoupled from the RG flow due to their finite mass. Therefore the 
influence of the gauge-sector on the (chiral) order-parameter potential (i. e. matter sector) 
in the deep IR is suppressed. Thus, chiral symmetry breaking is mostly sensitive to the 
running of the coupling in the mid-momentum regime rather than its running in the deep 
IR. In this respect, the findings in the present paper are in accordance with the findings in 
where the (de-) confinement phase transition in pure SU(2)- and SU (3)- Yang-Mills 



Ref. 



theories has been studied. 

As discussed in Sec. IIIDl we have also checked that the results are only sensitive to 
our ansatz for the momentum dependence of the four-fermion interaction at the percent 
level. Apart from these errors, systematic errors enter our calculations mainly from omitting 
instanton effects and deconfinement dynamics as described by the Polyakov-Loop. We expect 
that the transition temperature becomes larger when we include these effects. Thus the 
result for Tc given in Eq. (!73|) should be considered as a lower bound for the transition 
temperature. We address these issues further after the discussion of our results for the 
chiral phase boundary. 

In order to determine the curvature of the chiral phase boundary of 1-flavor QCD at small 



12 This can be seen from looking at the ratios k^,^{T = 0,// = Q)/k^f-{T = 0, = 0) « 0.70 and 
T^^{fi = 0)/T^'^^-{iJ. = 0) w 0.69. The almost perfect agreement of the two ratios indicates that the 
finite-temperature flows are less sensitive to the details of the running of the gauge coupling in the deep 
IR. 
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Author (s) Ref. 


Method 




Nf 

J 




J. Braun this work 


QCD RG flows 




1 


0.97.. 1.28 


J. Braun this work 


QCD RG flows: vs. Ocr 





1 


0.40. .0.47 


P. de For cr and et al. [2] 


Lattice QCD: imag. /i 


0.032 


2 


0.500(54) 


F. Karsch et al. |Y| 


Lattice QCD: Taylor + Rew. 


0.005 


3 


1.13(45) 


P. de For cr and et al. [5] 


Lattice QCD: imag. /i 


0.026 


3 


0.667(6) 



Table I: Comparison of the results for the curvature of the QCD phase boundary for A'^c = 3 as 
obtained from Lattice QCD and QCD RG flows. The table is not exhaustive and we concentrate 
here on results for degenerate quark flavors with (current) mass arric, where a denotes the lattice 
spacing. The RG calculations have been performed in the chiral limit. Apart from the current 
quark mass, the Lattice QCD simulations differ in the Lattice volumes and spacings used as well 
as in the implementation of the fermions. 

quark chemical potential, we compute Tc(/i)/Tc(0) for < fi/Tc{0) < 0.7. We then extract 
the curvature of the phase boundary at vanishing chemical potential from these results. At 
small chemical potential the phase boundary can be expanded in powers of /x^, which yields 

|^^l-M%^.mj(^)%^... (74) 

The coefficient depends on the number of quark flavors Nf, the number of colors Nc and 
the current quark mass. These dependences can be understood qualitatively by looking at 
the underlying mechanisms for chiral symmetry breaking as we discussed them in Sec. IIII Al 
The phase transition temperature and the curvature itself are sensitive to the magnitude of 
the quark condensate at T = and therefore also to the magnitude of the constituent quark 
mass at T = 0. The larger the constituent quark mass is for a given chemical potential /i, 
the smaller is the impact of fi on the phase transition temperature. Therefore becomes 
smaller with increasing (current) quark mass m^. 

The scale for the dynamically generated quark mass is essentially set by the scale kcr at 
which the gauge coupling exceeds its critical value. The scale fccr is related to the scale 
Aqcd where the gauge coupling strongly increases, thus kcr ~ Aqcd- The flavor- and 
color-dependence of Aqcd can be estimated from the position of the Landau pole of the 
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Figure 3: Phase boundary of 1-flavor QCD in terms of the dimensionless quantities Tcr{f^)/Tcr{0) 
and /i^/(7rTcr.(0))^. The hght gray band represents the results from the RG flows including the 
effects of gauge degrees of freedom. The black band indicates the result from an estimate of the 
phase boundary from a study of the running of the strong coupling as versus Ocr along the lines of 
the study in Ref. [31 



32( 1 . The width of the bands gives the theoretical error due to uncertainties 



in the gauge sector. 



perturbative one- loop running of the coupling^^ 32j: 

Aqcd ~ Mz e"3«M^ ~ Mz e~^^^^^^^ (l - eNj + C((eiV/)^)) (75) 

where 

We observe that Aqcd decreases linearly with Nf. Therefore the constituent quark mass 
becomes smaller with increasing Nf and bigger. For fixed Nj but increasing Nc, we find 
that Aqcd increases because In ^^^^ ~ Thus the quark mass increases with Nc and 
therefore decreases with increasing N^. A quantitative study of the Nc dependence of 



We neglect terms arising from the chemical potential, since Aqcd ^ A* foi' the values of the chemical 
potential we have in mind. 
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the phase boundary is work in progress [94] • We stress that these simple considerations of 
the dependences of the curvature in terms of the running coupling are in agreement with 
large A^c-considerations of the QCD phase boundary. In Ref. 95|, it has been shown that 
K,fj_ ~ Nf/Nc in leading order of an expansion in powers of l/Nc- 

In order to compare our results for the curvature of the phase boundary to results from 
lattice QCD simulations, we present our result for the phase boundary in terms of the 
curvature k^. In the first line of Tab. [T] and Fig. [3], we present our results for the phase 
boundary for QCD with one quark flavor together with results for two- and three-flavor 
QCD as obtained from lattice QCD simulations. We obtain 

n^^^- = 0.97 and /t™ = 1.28 (77) 

for ORef. and apE, respectively. From this, we conclude that the curvature depends signifi- 
cantly on the gauge-field dynamics whereas it is expected that the critical dynamics at the 
phase transition can be described with a simple 0(2) model due to universality. 

We list different lattice results for in Tab. [H From this table, we read off that the 
lattice result for for three degenerate quark flavors obtained with a Taylor expansion 
of the path integral is approximately twice as large as the result obtained with imaginary 
chemical potential. From our discussion above, we expect that decreases roughly linearly 
with decreasing Nf. Lattice QCD simulations with imaginary chemical potential are in 
agreement with this expectation. We compare our result for Kf^{Nf = l) with a naive linear 
extrapolation of n^{Nf = 2) and K^{Nf = 3) obtained from Lattice simulations with imaginary 
quark chemical potential^^ to A'^^ = L We then find that our result is roughly twice as large 
as the extrapolated value from these Lattice simulations. 

Let us finally discuss the systematic errors in our results arising from neglected operators 
associated with instantons and quark confinement. Although our truncation allows us to 
study chiral symmetry breaking with its underlying mechanisms in terms of quarks and 
gluons, it does not yet allow us to study the deconfinement phase transition. Within the 



RG 
2 



33 



ramework, the deconfinement phase transition in pure Yang-Mills theory with A'c 



3J], and with A'c = 3 [33|, has recently been successfully studied. The next step would 



Suitable results for obtained with sufficiently small quark masses for Nf = 2 are presently not available 
from Lattice QCD simulations in which a Taylor expansion combined with reweighting techniques has been 
applied. 
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be to couple the RG flow of the Polyakov to the present approach. We can estimate the effects 
of the inclusion of the Polyakov-loop by considering Landau-DeWitt gauge, D^{A){A— A) = 
0, where we identify the background field with the Polyakov-loop field, A = (Ao). The quark 
propagator effectively acquires an additional mass term since the vacuum expectation value 
of the zero-component (All of the gauge-field shows up in the Matsubara frequencies of the 



quark fields as follows^^ 20[ |: 



K + i/i + ^(Ao))^ (78) 

As we have discussed above, a larger quark mass translates into a flatter curvature (smaller 
K^) and a higher (chiral) phase transition temperature, so this would be the likely effect of 
the Polyakov loop as well. An inclusion of the Polyakov-Loop dynamics in our RG study 
would eventually allow for a dynamical study of deconfinement and chiral phase transition 
and their interplay at the same time. 

Instanton effects associated with the Ua{^) anomaly are certainly important for a more 
quantitative prediction of the curvature of 1-flavor QCD phase boundary since they directly 
influence the quark propagators. Our results for the phase boundary have been obtained 
from an ansatz for the effective action which has a global Ua{^) symmetry. This symmetry 
is broken in QCD by the presence of gauge-field configurations with non-trivial topology. 
Instantons are such gauge-field configurations. In the context of instantons, the ?7 ^(1) 
symmetry is broken by their induction of masslike fermion interactions in 1-flavor QCD^^ 
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47 



48|: 



451, 



d^xmj{tljRtljL - ^l^r) = j d^xmj (^/'Ts^/') • (79) 

The associated mass parameter mj is exponentially suppressed for small gauge coupling 
Cisi fT^T I Anqn ~ e~^'^/"% but becomes significantly large near the chiral transition scale, see 



e.g. Ref. [30|. Due to the presence of such an instanton mediated interaction, the pions ac- 
quire a mass in the deep IR as well. Therefore the instanton-mediated interactions influence 
directly the curvature of the phase boundary and the phase transition temperature. Since 
the instanton-mediated interactions act like an explicit mass term for the quark fields, we 



15 



We have (^o) = only for T — > oo. 



The treatment of instantons within the functional RG framework has been discussed in detail in Refs. 
96|. 



3C 



35 



expect the curvature to become flatter, the phase transition temperature to become higher 
and the second order phase transition to turn into a crossover. 

Although we do not study the effects of a broken global Ua{^) symmetry and an inclusion 
of the Polyakov-Loop dynamics explicitly, we can give an estimate how they affect the 
phase boundary by using Eq. ( 1701) to determine the phase boundary. Estimating the phase 
boundary from Eq. (1701) means to effectively stop the RG flow at the scale kcr at which the 
strong coupling exceeds its critical value. We flnd k^r ~ 1.5 GeV for ttRef. and kcr ~ 0.9 GeV 
for apE- Our results for the curvature are given in Tab. [B We flnd that the estimated 
values for are 



K 



J"^- = 0.47 and = 0.40 , (80) 



which are smaller by roughly a factor of two compared to those values obtained from a solu- 
tion of the full set of flow equations, see Eq. (1771) and Tab. HI This is due to the fact that the 
massless Goldstone modes dominating the IR physics are effectively cut off when we estimate 
the phase boundary from Eq. fl70l) . Instanton- mediated fermionic interactions effectively in- 
troduce an IR cutoff for both the fermions and the bosons in the IR which is roughly of the 



order of 1 GeV in case of 1-flavor QCD 30|. Therefore our estimate of the curvature of the 
phase boundary from Eq. ( 170ll is less contaminated by these topological aspects of QCD and 
serve as an estimate for a lower bound for the curvature k^. A quantitative analysis of the 



influence of topological effects on the QCD phase boundary is postponed to future work [97|. 
Moreover we observe that the uncertainty for the curvature in Eq. (IHOj) . obtained from our 
estimate for the critical gauge coupling, is much smaller than the uncertainty found in the 
results from a solution of the full set of RG flow equations, see Eq. (|77j) . This is due to 
the fact that the running coupling exceeds its critical value on scales where its running is 
still close to the perturbative 2-loop running and the differences between aRcf. and apE are 
therefore small. On the other hand the results (177|) for the curvature from the solution of 
the full set of RG flow equations are sensitive to the non-perturbative running of the gauge 
coupling in the mid-momentum regime {p ~ 0.5 — 1 GeV). Therefore the uncertainty in the 
curvature reflects mostly the uncertainty in our truncation of the gauge sector. 
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IV. SUMMARY AND CONCLUSIONS 

In this paper we have presented a functional RG approach which allows to study the 
QCD phase boundary from first principles. Our work aims to set the stage for future studies 
in this direction incorporating two and three quark flavors. As a first application, we have 
computed the phase boundary of 1-flavor QCD at small chemical potential and found a 
second order phase transition if Ua{^) violating terms are neglected. However, our study 
in its present form does not allow us to detect a critical endpoint in the phase diagram of 
1-fiavor QCD since it relies on a low-order expansion in n-point functions in the scalar sector 
with n < 4. Therefore the reliability (i. e. the radius of convergence) of a Taylor expansion 
of the chiral phase boundary in powers of the quark chemical potential cannot be checked. 
With respect to the nature of the phase transition, however, we expect our truncation to 
be reliable for small chemical potentials (at least for // = 0). In this regime the nature 
of the finite-temperature phase transition is dominated by the underlying 0(2) symmetry 
while quark effects are subleading^''. An improvement of our truncation with respect to an 
inclusion of higher n-point functions in the scalar sector of our truncation, which allows us 
to search for a critical endpoint, is deferred to future work. 

Apart from a numerical study of the phase boundary, we have discussed the underlying 
mechanisms of chiral symmetry breaking in terms of quark-gluon dynamics and how these 
mechanisms relate to the dependence of the curvature on the number of quark flavors A'^^, 
the number of colors Nc and the current quark mass rric. In particular, we have argued that 
the curvature is linearly dependent on Nf/Nc- 

Although the present approach contains already all ingredients that are necessary for a 
study of the QCD phase diagram from first principles, our numerical predictions for the 
phase boundary at small chemical potentials suffer from the underlying approximations. 
As we have already discussed above, our present low-order expansion of the scalar sector 
in terms of n-point functions does not allow us to detect a first-order phase transition. 
Moreover, we have further theoretical errors entering our study due to our truncation of 
the gauge sector. In order to "measure" the uncertainties arising from this sector we have 
used the running gauge coupling from a functional RG study employing the background- 



Even though gluon-induccd quark interactions drive the system towards the phase transition, quark effects 
can be suppressed at the pliase transition. 
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field method 



31 



32l | and from lattice QCD 77|] as an input. Thereby we have exploited the 



fact that the coupling of the gauge sector and the matter sector is dominantly given by the 
running gauge coupling, the wave-function renormalizations of the gluons as well as the quark 
mass. We found that the curvature of the phase boundaery is about 30% smaller for the 
lattice coupling than for the background-field coupling. This suggests that the curvature of 
the phase boundary is indeed sensitive to the underlying gauge-field dynamics beyond large 
Nc which is an important result for presently used (P)NJL-type models. With respect to the 
current debate on scaling versus decoupling scenario in the gauge sector, see e. g. Ref . [93| , 
our analysis suggests that the behavior of the propagators in the deep IR does not strongly 
influence the shape of the phase boundary. This observation is accordance with a study of 



the deconfinement phase transition in pure Yang-Mills theory 2 



Apart from our analysis of the role of the gauge-field dynamics essentially stemming from 
the gluonic two-point function, we have an uncertainty originating from the fact that we 
have expanded the gauge-sector about a vanishing (Aq) instead of taking its temperature 



and scale-dependence into account |33|]. Concerning f/A(l)-violating terms, we have provided 
an estimated of their influence on the phase boundary. We have argued that we expect the 
curvature of the QCD phase diagram at small chemical potentials to become flatter by a 
factor of two when these two missing pieces are included. This lead us to an estimate for 
a lower bound for the curvature of the phase boundary. We add that this estimate for the 



curvature compares nicely with a linear extrapolation to Nj = 1 of lattice QCD studies for 
Nf = 2,3,4 with imaginary chemical potential, see Refs. fly, 3,3. Finally we would like 



to mention that instanton effects as well as an expansion of the gauge sector about a scale 
and temperature-dependent (An) can be included dynamically in our study by combining it 



with the findings in Refs. 



96| and |33|, respectively. 



We conclude that our approach allows us to compute the phase boundary unambiguously 
in the sense that the scale for all of our results is set by a single input parameter, namely the 
value of the strong coupling as at e. g. the Z-boson mass scale. This is a great advantage 
compared to studies of the QCD phase diagram in terms of (P)NJL-type models, in which 
the results for the curvature and the location of a (possible) critical endpoint in the phase 
diagram strongly depend on the choice of the UV cutoff and a set of input parameters. 
Therefore we think that the present approach is very promising since it allows to bridge the 
gap between quarks and gluons on the one hand and hadronic degrees of freedom on the 
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other and opens up the possibihty of studies of the QCD phase diagram with two and three 
quark flavors which depend on only a single physical input parameter. 
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Appendix A: DEFINITION OF THE PROPAGATORS 

We define the boson propagator as follows: 

Psipo, {Pi}) = -1, ^ • (Al) 

For the fermions, it is convenient to define the modified four-momenta: 

The fermion propagator can then be written as 

Pi;\M^) = V^^\M^) - (1 + T^r'M^] , (A3) 
pI~\m^) = P(")(M^) [^(") + il+r^r'M^] (A4) 

with 

^ -Zi{TPo,{TPi}){l + r,,) 

[zliTPo, {TPi}) {TPo + if^W + (ZiiTPo, {TPi})p{l + r^W + M^' 
We should keep in mind that is in general a complicated function depending on the 
background fields, the Yukawa coupling, the gauge coupling as well as the four-momentum^^: 

= M^po, {pi}, {$,}, h, g, {A^}) , (A6) 
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Note that the eouphngs depend in general on the momenta as well. 
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where $ represents a two-dimensional vector of real scalar fields. 
Finally, the gauge-field propagator is given by 

p ( i .X) ^ ( ^ P^^P 

" 4(Po,te})p^ + ^^(j>o,te})]?(l+rB) V ~¥ 



where ^ is the gauge-fixing parameter. 



(A7) 



Appendix B: THRESHOLD FUNCTIONS 

The regulator dependence of the flow equations is controlled by (dimensionless) threshold 
functions which arise from Feynman graphs, incorporating fermionic and/or bosonic elds. 
Let us first introduce the so-called dimensionless regulator-shape function rsix) and r^(x) 
for bosonic and fermionic fields. These functions are implicitly defined by the regulator 
function Ri as follows: 

Rb{Po, {P^}) = 4fc(Po, {p^})p'rB{pVe) (Bl) 
for the bosonic fields and 

RAPo, {P^}) = Zi,{po, {P^mr4p'/k') (B2) 
for the fermionic fields. For the gauge fields, we choose 

Ra{Po, {p.}) = 4.(Po, {p.})p'rs{pVe) (5,. - ^) + -^p'Mp'/k') . (B3) 



In this work, we employ a 3d optimized regulator-shape function 23|, |5l|, l52] : 



rB(x) = 0-1^ e(l-x) and = (^-^ - 1^ 0(1 - x) . (B4) 

In the following, we use these regulator shape functions whenever we evaluate the integrals 
and sums in our general definitions of the threshold functions. 

In order to define the threshold functions, it is convenient to define dimensionless prop- 
agators for the scalars (B), gluons (A) and fermions respectively: 

GBA^o,i^) = ^ — ^ ^ — (B5) 
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and 

where zb = z\/Zj^, za = z\/Zp and z^ = Z^!^/Z^ give the ratio of the wave- function 
renormahzation in the direction parallel and perpendicular to the heat-bath. 

First, we define the threshold functions which encounter in the flow equations for the 
effective potential. For the bosonic loops, we find 



n=—oo 



' ' + (B7) 



d- 1 v^TTTJ V + ly V2 

where t = T/k denotes the dimensionless temperature and u = 2Tmt denotes the dimen- 
sionless bosonic Matsubara frequencies. The function ub represents the Bose-Einstein dis- 
tribution function 

Bosonic threshold functions of order n are then derived from Eq. flBTD by induction: 

^/i^)'^'^) (t, t^b) = -{n + Sn,o) iSi^"^ {i, co; Vb) • (B9) 
For the fermion loops contributing to the flow equations of the effective potential, we find 

n=-oo 

" d- 1 y/TT^ ~ "rf^) ~^^(^~'/^'^) ~^V'(^~^-/^^^)) • (BIO) 

Here, we have introduced the dimensionless fermionic Matsubara frequencies z/„ = (2n+l)7rt, 
the dimensionless chemical potential fi = fi/k and the Fermi-Dirac distribution functions n^: 



1 



n4ir^,u;) = ^ eC^-VlT^). (Bll) 

Higher-order fermionic threshold functions are also given by induction: 

^^l^P'^'Kl ^, /i; m) = -{n + 5„,o) l^n+i'\l u;, fi- v^) . (B12) 
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Let us now define the threshold functions which are involved in the computation of the 
Yukawa coupling. We have 

L[^f^''''^\i, uj^, /i, ujb, Qo; Vi^^ Vb) 

= -- ^ / dxx'^dtG^{{Dn + ifif,uj^)GB{ii>n-Qoy,ujB). (B13) 

n=— OO 

In order to evaluate the integral over x (spatial momenta), we use^^ 

9, 



9t 



(B16) 



where the first and the second line tells us how the formal derivative dt acts on fermions and 
bosons, respectively. The threshold function reads then 

-^1,1 ''''' \ii^ipil^i^BiQQ'i'ni}iVB) 

'^'^ ^GB{{i>n-Qof,u;^)}, (B17) 



where we have introduced the auxiliary functions 

g^{xo,u}) = — , (B18) 

Z^Xo + 1 +UJ 

gB{xo,uj) = ^7^. (B19) 

ZbXq + 1 + LJ 

In Landau gauge (or any other gauge with gauge-fixing parameter ^ 7^ 1), we encounter 
additional threshold functions in the flow equations for the Yukawa coupling due to the 



Here, we give only explicit expressions for the formal derivatives obtained with the regulator shape func- 
tions (lB4l). 
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presence of longitudinal gauge bosons: 



1^^, /i, (^B, Qo'i VipJ Vb) 

i 

~ ~2 



~ oo ^ d-l 

/ dx 



n=—oo 
oo 



[i>n - QgY + X 



n=— oo 



+(^{.%((^n-gof) 



where the auxiliary functions Hd{z) are given by' 



20 



1 / 

dx 

\z + x 



(B20) 



(B21) 



Next, we discuss the threshold functions which are involved in the computation of the 
wave-function renormalizations. We shall start with the threshold functions for the fermionic 
wave-function renormalization due to its close relation to those for the Yukawa coupling. 
Overall, we have three different types of threshold-functions contributing to the flow of Z^. 



The first one is closely related to those functions found in Refs. 

.{FB),(d)ir 



30 



62| and reads: 



M\2 \t,UJ^,^l,UJB,Qo]y]^,TlB) 

~ oo d 1 ~ ( ~ d ~ ~ ~1 

= 2 ^ / ^^^^^t{(l + '^V')G'v>((i^n + i/i)',c^^)^G'B((z>„-go)',^B)j-(B22) 

n=— oo 

Evaluating the integration over x, we find 

j^{FB),{d)^^^ ^V" ^B, Qo) 

oo 2 

n=— oo 

Due to our choice of the regulator functions, is independent of Z^, even for 

dtZ^ 7^ 0. The second type of threshold function that we encounter is given by 

,uj^, fi,ujB,Qo','r]^,VB) (B24) 

^ ^ roo ~ ( ~ d ~ ~ ^ 

= o 511 / ^^7: 5^Ul + r^)Gv,((^>n + i/i)^t^V■):r'^^((^"~^o)^^^B) 

^ ^t^^ Jo [vn - QoY + X I- dx 



20 



The functions Hd{z) are related to the Hypcrgeomctric function ^Fi^ see e. g. [98| 
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As in the case of the threshold functions for the flow of the Yukawa-coupling, this threshold 
function as well as the next threshold function that we are going to discuss is only present 
in gauges with gauge- fixing parameter ,^ 7^ 1. Due to the similar momentum structure of 
Ji4^^^^'^'^\ we immediately obtain 

= i Y] Gi>{{i>n + i/i)^^^) (GBiii^n - (5o)^^B)) ( — — 

- Y^yH(^n - Qo)^)) (B25) 
for the given choice of the regulator functions. The last threshold functions present in the 



RG flows of Zjjj is a generalization of an already know threshold function found in Ref. 



30|: 



, uj^, yU, lub, Qo', Vb) 



(B26) 



4E 



dx- 



d-3 
X 2 



2 ^t^^Jo {{un-Qoy+xy 



dtUl + r^)G'^((z/„ + i/i) , uj^)GB{ii^n - Qo) , ^^b) 



The integral over x can be performed straightforwardly, yielding 

, uj^, /i, ujb, Qo', Vi>, Vb) 



n=— 00 

-V^ [<%((^n-go)^) - HlX!,^{{u^-Qo, 



2^^ ( (z>„ + i/i) ^ t^^)i^(24) ( (i^„ - Qo)^) 



+ 2ri^g4ii>n + r^)^C0^,) terJ4)((^"-^0)') - H[^^3)ii^n-Qo)') 



Am) 



-H, 



{T-2)(('^n-Qo)^) 



(B27) 



The threshold functions needed for the computation of the scalar anomalous dimensions 
are given by 

i °° 



Mi^^'^'^\i, uJi,uj2;r]B) 



f; l^^dxx'^dt (^-^GB{Colu,)^(^-^GB{Ci;luj,)^ .(B28 



The integration over x can be performed analytically for the regulator shape functions under 
consideration and we obtain 



(B29) 
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Note that is independent of 77^ for our choice of the regulator functions as it has 

already been found in Ref . 99|] . Two more contributions to the scalar anomalous dimensions 
come from purely fermionic loops. One of those contributions is proportional to the square 
of the vacuum expectation value of the scalar field and the corresponding threshold function 
reads 

2 



-- ^ y dxx'-^dt\^—G^{{un + ijj)\uj)^ . (B30) 



The integration over x can be carried out analytically and we find again that ' 
independent of the fermionic anomalous dimensions for the cutoff function under consider- 
ation: 



(B31) 



The threshold function of the second contribution consisting solely of fermionic internal lines 
is given by 



M 



{F),{d),7 



t 



dxx 2 dt 



d 

dx 



(l + r^)G^((^>„ + i/i)^cu) . (B32) 



Performing the integration over x, we find 



(B33) 



2d-6 4^ 

Finally, we have to discuss the threshold functions corresponding to 1 PI box diagrams. 
These threshold functions are needed for a computation of the RG flow of the four-fermion 
interaction. As in the case of the Yukawa-coupling and the fermionic wave-function renormal- 
ization, we encounter additional threshold-functions in gauges with gauge-fixing parameter 
^ 7^ 1. We start, however, with the discussion of the threshold functions that are closely 



related to those already discussed in Refs. (29|, l30|, Il00( | : 



j(FB),{d)i7 ~ ~ \ 

-^l,l,ni,n2l''' ^i'l A'l' /^2, ^B,!-, ^B,2-, Vi>^ Vb) 



■ E 



dxx 2 



X GB(i>l,UJBA 



X 



ni 



GB{yl,u;B,2) 



n2 • 



. (B34) 
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By performing the integration over x, we obtain the expression for the threshold function 
used in this work: 

- iFB),(d) 



1,1, rti ,n2 

(t, w^, /ii, /i2, ujb,i,ujb,2; V^, Vb) 



2t 



d 



_ — -\ r 1 "1 r 



"2 



X 



Since we are exphcitly studying Landau-gauge QCD, we have an additional class of threshold 
functions given by 

£iS'""(*.'^*.Ai,A2.^B;%,^B) 



■ E 



d+1 



X '2 

dx— dt 

u^+x 



(l+r^)G^((^/„ + i/il)^Cl;^,) x 
(l+r^)G'^((^>„ + i/i2)^u;^) GBivl^us) }> . (B36) 



Inserting the regulator shape functions (lB4p . we find 

4S'^^^(t,-.,/ii,/i2,-B;./„r/.) 

oo 



X 



'(d-i)l 
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